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Abstract 


Open rotors are once again under consideration for propulsion of the future airliners because of 
their high efficiency. The noise generated by these propulsion systems must meet the stringent 
noise standards of today to reduce community impact. In this paper we review the open rotor 
noise prediction methods available at NASA Langley. We discuss three codes called ASSPIN 
(Advanced Subsonic-Supersonic Propeller Induced Noise), FW - H pds (Ffowcs Williams-Hawk- 
ings with penetrable data surface) and the FSC (Fast Scattering Code). The first two codes are in 
the time domain and the third code is a frequency domain code. The capabilities of these codes 
and the input data requirements as well as the output data are presented. Plans for further 
improvements of these codes are discussed. In particular, a method based on equivalent sources 
is outlined to get rid of spurious signals in the FW - H pds code. 
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1- Introduction 

The high fuel prices of recent years have caused the operating cost of the airlines to soar. In an 
effort to bring down the fuel consumption, the major aircraft engine manufacturers such as GE, 
Rolls-Royce and CFM International are now taking a fresh look at open rotors for the propulsion 
of future airliners. Open rotors, also known as propfans or unducted fans, can offer up to 30 per 
cent improvement in efficiency compared to high bypass engines of 1980 vintage currently in 
use in most civilian aircraft. The lower fuel consumption of open rotors also contributes to the 
reduction of carbon emission into the atmosphere by airliners. Open rotor propulsion system can, 
therefore, be considered a Green engine. The single and contra-rotating open rotors are both 
under consideration for powering future aircraft. An example of a contra-rotating open rotor is 
shown in Figure 1 below. 



Figure 1- An example of contra-rotating open rotor design by GE 

Open rotors operate at transonic to supersonic helical tip speeds that can go as high as Mach 1.3. 
The noise of these propulsion systems can cause both passenger and community concerns. The 
ability to predict the noise of open rotors at the design stage and to reduce the noise to an accept- 
able level are, therefore, of great importance to engine manufacturers. In 1980s, a comprehen- 
sive study of open rotors by NASA and the U.S. aircraft engine industry was carried out involv- 
ing theory development, model and full scale tests. An impressive array of technology develop- 
ment in aeroelasticity, composite materials, mechanical design, aerodynamics and aeroacoustics 
resulted from this study. A large number of technical reports and publications in open literature 
were written by the engineers and researchers involved in this effort. These reports and publica- 
tions are of great value to those involved in the resurgence of interest in open rotors around the 
world. 
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NASA Langley researchers have contributed significantly to the development of aeroacoustic 
technology of open rotors. The present authors feel that it is appropriate at this time to present 
the current noise prediction technology at Langley and review the input data requirements, 
strengths and limitations of each method as well as the associated problems in need of attention 
by the researchers. These are the aims of the present paper. 

In Section 2, we present a brief history of research on the aeroacoustics of rotating blade machin- 
ery at Langley Research Center. In Section 3, we discuss the available noise prediction codes for 
open rotors developed at NASA Langley and their capabilities. In particular, we briefly talk 
about the two useful formulations used for the computation of noise from subsonic and super- 
sonic surfaces. Here we discuss the open rotor noise prediction codes ASSPIN and one based on 
Ffowcs Williams-Hawkings equation with penetrable data surface (FW - H pds ). The scattering 
of sound from surfaces near the rotor are calculated using the fast scattering code (FSC) which is 
also discussed in this section. Plans for further improvements of these codes are given in Section 
4. Here we propose a method based on equivalent sources to establish the null field inside the 
data surface to get rid of the spurious signals in the FW - H pds code. 

2- A Brief History of Research on the Aeroacoustics of Rotating Blade 

Machinery at Langley 

Langley has been in the forefront of research on aircraft noise generation and prediction since 
the 1930's [SI]. Some of the early researchers in this subject at Langley were Deming, Garrick, 
Watkins, Hubbard and Maglieri. This research was aided greatly with the introduction of the 
high speed digital computers in the 1960's that were used both in acoustic experiments as well as 
in the development of noise prediction codes. About the same time, the pace of research in the 
field of aeroacoustics in universities, industry and government agencies accelerated due to the 
public objection to aircraft noise both in the cabin and around the airports. Modern digital com- 
puters have made it possible to use advanced mathematical models to study the acoustics of 
propulsion systems with realistic geometry and operating conditions. NASA scientists and engi- 
neers have contributed significantly to aircraft noise reduction both by theory development and 
the use experimental methods. 

In what follows, we will discuss the history of theoretical model development at Langley that the 
author was involved in the period 1975 to present. We will confine ourselves to subjects of 
relevance to aircraft with open rotors as their propulsion system. This would also involve the 
calculation of installation effects and the noise foot prints at various flight conditions. 

The introduction of the concept of the acoustic analogy (AA) by Lighthill in 1952 [S2] was one 
of the most important advances in aeroacoustics which has resulted in current high level of 
aircraft noise prediction technology. In 1969, Ffowcs Williams and Hawkings presented an 
extension of Lighthill's jet noise equation to the noise from turbulence in the presence of moving 
surfaces appropriate for studying the acoustics of rotating machinery- helicopter rotors, 
propellers, fans and open rotors [S3]. This extension is known as the Ffowcs Williams- 
Hawkings (FW-H) equation. 

Let a moving surface be described implicitly by fix, t) - 0 such that Vf = n, n is the unit out- 
ward normal. This assumption implies that / > 0 outside the moving surface (see Figure 2). 
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f(x,t)< 0 


Figure 2- The definition of the moving surface implicitly as f(x, t) - 0. Note that v / = n 
where n is the unit outward normal to the surface. 


The Ffowcs Williams- Hawkings (FW-H) equation is 
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Here and elsewhere in the paper, we will use the summation convention on the repeated index. 
In this equation D 2 is the wave operator in three dimensional space. Also p' = c 2 p' = c 2 (p - p 0 ), 
c and po are the speed of sound and density in the undisturbed medium, respectively. Note that 
p' can only be interpreted as the acoustic pressure if p'/Po« 1- The symbols v„, p and 
Tjj - p Ui iij - CTjj + (p' - c 2 p') djj are the local normal velocity of the surface, the local gage 
pressure on the surface (in fact p - p 0 ), and the Lighthill stress tensor, respectively. In the defini- 
tion of the Lighthill stress tensor, cry is the viscous stress tensor and 6y is the Kronecker delta. 
The Heaviside and the Dirac delta functions are denoted H(f) and 6(f), respectively. In the 
second term on the right of eq. (1), we have neglected the viscous shear force over the blade 
surface acting on the fluid. 

The three inhomogeneous source terms on the right of Eq. (1) are known as the thickness, load- 
ing and quadrupole sources, respectively. Note that the thickness and loading sources are sur- 
face sources while the quadrupole source is regarded as a volume source although it can be 
analytically broken into volume and surface sources [21, 23, 27, 31, 36]. This latter point of 
view can explain several (old and new) mechanisms of noise generation in rotating blades some 
of which has not been studied in detail, e.g., flapping wakes in the vicinity of trailing edge. 
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We note that we have artificially converted a nonlinear problem of noise generation by a moving 
surface to a linear problem by using the acoustic analogy. All the nonlinearities are lumped into 
the Lighthill stress tensor which is assumed known from near field aerodynamic calculations. 
When we started working on helicopter and propeller noise in the early seventies, because of the 
limitations of digital computers, the most we could expect from aerodynamic calculations was 
the blade surface pressure. For this reason, using some physical reasoning, we neglected the 
quadrupole volume sources in FW-H equation and concentrated on development of formulations 
for the prediction of thickness and loading noise. Later on, as computers became more powerful, 
we included quadrupoles in our noise prediction. There was, however, another theoretical 
advance which led to the use of purely surface sources to which we will turn next. 

It was Ffowcs Williams himself who proposed to use a penetrable (porous or permeable) data 
surface to account for nonlinearities in the vicinity of a moving surface. We again assume that 
the penetrable surface defined by f (x, t) - 0 and the fluid velocity is denoted by u. The FW-H 
equation for penetrable (permeable, porous) data surface, FW-H pds , is: 


(3 <3 ^2 
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We have used the following notations in the above equation: 
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Li - p 8 ij fij + p Uj (u n - v n ) ( 4 ) 

where djj is the Kronecker delta. As in the case of eq. (1), in the first term on the right of eq. (4), 
we have neglected the viscous shear force over the data surface acting on the fluid exterior to the 
surface. The philosophy behind using FW-H pds is to locate the data surface / = 0 to enclose a 

moving surface, in such a way that all quadrupoles producing non-negligible noise are included 
within this surface. Therefore, no volume integration of the quadrupoles outside the data surface 
is necessary. High resolution CFD calculation is performed in the near field region (including 
turbulence simulation if broadband noise prediction is required). The data surface used for the 
acoustic calculation should be located within the region of high resolution CFD computation. As 
discussed in Section 3, the optimal location of the data surface, particularly when vortices cross 
the surface, is still the subject of research. One would like this surface to be as small as possible, 
because of the computer intensive nature of aerodynamic and turbulence simulation that require 
fine grid sizes and small time steps. 

We have developed several equivalent solutions of the FW-H equation in the time domain 
(called formulations in this paper) at NASA Langley that are used to develop two efficient 
codes for the prediction of the noise of rotating machinery. In the next section we present a 
discussion of these formulations, computer codes based on them and their capabilities. Note that 
the Fast Scattering Code (FSC) does not depend on these formulations but is included in the next 
section because we will find it useful for open rotor noise prediction when scattering of the 
generated noise by fuselage and wing surfaces becomes important. 
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3- Noise Prediction Codes for Open Rotors Developed at NASA Langley and 

Their Capabilities 

The first author has derived several solutions of the FW-H equation that are useful for efficient 
numerical computation of the noise of rotating blades. We mention here three formulations in 
the time domain valid in the near and far fields that include only the thickness and loading 
sources of the FW-H equation. In chronological order, these are: 

1- Formulation 1 [1, 8, 52]: This was developed for surfaces moving at subsonic motion if used 
with the Doppler singularity. It is valid for surfaces moving at transonic and supersonic speeds if 
at no time during the emission process the condition M r - 1 occurs and the collapsing sphere 
form of the formulation is utilized [1], Here M r is the Mach number of the source in the radiation 
direction. Note that in this section, the Mach number is always based on the speed of sound in 
the undisturbed medium. In the case of transonic and supersonic speeds, the multiple emission 
times of the source must be accounted for in the acoustic calculation. Formulation 1 is the sim- 
plest of our three formulations distinguished by an observer time derivative before one of its two 
surface integrals. This time derivative must be evaluated numerically which can be a source of 
numerical error. 

2- Formulation 1A [7-9, S4, 52]: This result was obtained by taking the observer time deriva- 
tive of Formulation 1 inside the integral analytically and, thus, eliminating the need for taking 
the derivative numerically. Although, in principle, this result can be obtained by elementary 
calculus (not recommended), the simplest method of deriving it is using an identity of general- 
ized functions for the free space Green's function of the wave equation at the beginning of the 
derivation [8]. An immediate consequence of eliminating the observer time derivative of Formula 
tion 1 is the increase in the complexity of the formula and the increase in the order of the 

Doppler singularity from (1 - M r )~ l to (1 - M,.) -3 . The increase in the complexity of the formula- 
tion does not cause any concern for numerical implementation on the computer. But the increase 
in the order of the singularity increases the sensitivity of the result near M r - 1 . However, in 
practice, we must switch to Formulation 3 in the acoustic code when this situation occurs. The 
conditions for the validity of Formulation 1A are similar to Formulation 1 and is used in the 
forms with Doppler singularity (subsonic sources) as well as collapsing sphere (transonic and 
supersonic sources). An important use of Formulation 1 A is in the acoustic code based on the 
FW - H pds equation. 

Formulation 3 [13, 16, 17, 20]: This complicated result is very difficult to derive and requires 
advanced knowledge of differential geometry, general tensor analysis and generalized functions. 
The formulation is valid at all range of the speed of the source surface- subsonic, transonic and 
supersonic. However, because of its complexity, it is recommended to use it for the parts of the 
surface for which the condition M r - 1 occurs. For blades with sharp leading and trailing edges 
in the supersonic region of the blade and properly designed sweep, no singularities should 
appear in the acoustic pressure when this formulation is used. Otherwise, one can expect the 
integrable logarithmic singularity in the acoustic pressure [22, 27], 


FF, MHD, AT& DMN- AIAA 2009.nb | 7 


An explanation of why Formulation 3 is so complicated is given in [46]. An important test of 
this formulation is to compare the predicted acoustic pressure signatures and spectra from a 
subsonic rotating blade exclusively from Formulations 1A and 3 under identical input data. This 
computation has been performed at NASA Langley, documented but not published. The agree- 
ment of the predicted results from the two formulations was excellent for both the thickness and 
loading noise. Furthermore, this agreement also verified that both these formulations were coded 
correctly in ASSPIN. 

Some old and obsolete NASA Langley acoustic codes used Formulation 1 for high speed pro- 
pellers [6]. We currently use only Formulations 1A and 3 in our noise prediction codes such as 
ASSPIN discussed in the next section. The best article to read about how to implement these 
formulations most efficiently to predict the noise of rotating blade machinery is [20]. 

We now turn to the description of two computer codes for open rotor noise prediction. Note that 
our codes are valid for an observer in the near and far fields. The third code we discuss below 
is the Fast Scattering Code (FSC) which is recommended for the prediction of the effect of the 
scattering of the open rotor noise by fuselage and wings of an aircraft. This phenomenon 
includes the important shielding effect of noise by the rigid surfaces of the aircraft known also as 
the installation effect. 

3.1- Advanced Subsonic-Supersonic Propeller Induced Noise (ASSPIN) Code 

This code was developed in mid-nineteen eighties originally for a single row of open rotors 
known as propfans that operated at supersonic tip speed at cruise condition [13, 16, 17, 20]. It 
can also calculate the noise of a contra-rotating propeller with equal number of blades in each 
row. However, the present code can form the core of a much more capable open rotor noise 
prediction code by the addition of pre- and post-processors. There is an excellent Users Manual 
written by Dunn and Tarkenton [S5]. This code uses Formulation 1A and 3 in such a way that 
the latter formulation is used for panels on the blades with near sonic speed. Farassat, Padula and 
Dunn [20] have described the structure and the algorithms of ASSPIN in great detail with some 
comparison of the predicted and measured acoustic data for propfans. Other publications of 
Langley with more examples of application of the code are [30, 34, 35]. One of the most signifi- 
cant tests of ASSPIN was against the frequency domain code of Donald B. Hanson with identi- 
cal blade geometry, operating conditions and blade aerodynamics. The results with good agree- 
ment were published in the same issue of the AIAA Journal in 1992 [35, S6] 

The ASSPIN code was written in FORTRAN 77 and developed originally on a CDC-7600 main- 
frame computer of NASA Langley. It runs very efficiently on a modem scientific desktop com- 
puter. Noise footprint computations can be performed routinely in a matter of few hours and 
even the animation of the footprint during the takeoff condition can be simulated in a reasonable 
time. Other than the fact that the nonlinear effects are not included in this code, it is one of the 
most robust, accurate and verified codes in the time domain in U.S. and, to our best knowledge, 
in the world. 

3.1.1- Input Data to ASSPIN 

1- Hot blade geometry- several options are available to users 

2- Number of blades in each row 

3- Inflow angle to the rotor axis 

4- Unsteady blade surface pressure 
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3.1.2- Output Data From ASSPIN 

1- Acoustic pressure time history for a stationery observer with respect to the ground or moving 
with the aircraft at uniform speed 

2- Acoustic pressure spectrum for an observer moving with the aircraft at uniform speed or the 
acoustic spectrum for contiguous segments of time for a stationary observer (nonperiodic time 
history). 

3.2- Ffowcs Williams-Hawkings with Penetrable Data Surface (FW - H pds ) Code 

It has been known since the early days of high speed propeller noise research that nonlinear 
propagation effects in the near field are important. To include the nonlinearities, as discussed in 
Section 1, the use of a prediction code based on the FW - H pds is an obvious choice. Since all 
the commercial airliners will fly at subsonic cruise speed in the foreseeable future, Formulation 
1 A is the appropriate solution of the FW - H pds equation to use. For an aircraft in steady rectilin- 
ear flight (the most common situation of interest), there is a simple closed-form formula based 
on the Garrick triangle construction for the computation of the retarded time in the acoustic 
calculations. Thus, the development of a noise prediction code based on FW - H pds equation 
using Formulation 1 A for this case is a relatively simple task. Nark and Farassat have developed 
such a code at Langley and have used it for ducted fan noise prediction [50]. It can be used for 
open rotors with some additional pre- and post-processors. 

When turbulence and wakes from rotors cross the data surface, spurious acoustic signals appear 
in the acoustic signature which reduce the accuracy of this method. We will discuss in Section 3 
our plans to alleviate this shortcoming of the method. 

3.2.1- Input Data to FW - H pds Code 

High resolution unsteady CFD data (static pressure, fluid density and velocity) on the data sur- 
face must be specified. 

3.2.2- Output Data From FW - H pds Code 

1- Acoustic pressure time history for a stationery observer with respect to the ground or moving 
with the aircraft at uniform speed 

2- Acoustic pressure spectrum for an observer moving with the aircraft at uniform speed or the 
acoustic spectrum for contiguous segments of time for a stationary observer (nonperiodic time 
history). 

3.3- The Fast Scattering Code (FSC) 

This code was developed at NASA Langley by Dunn and Tinetti based on the equivalent 
source method (ESM) [S7, S8]. The theory behind the FSC is explained well in these refer- 
ences. It is based on the equations of time harmonic, linearized acoustics and the ESM is used 
for solving an exterior Helmholtz equation as a boundary value (BV) problem for a surface 
moving at uniform rectilinear speed. The incident sound source, e.g., the open rotor noise, is 
Fourier analyzed and the FSC is then used to compute the scattered noise for each frequency 
component. The ESM is a much more efficient method of finding the scattered noise than other 
well known techniques such as the boundary element method (BEM). The advantage comes 
from the fact that in ESM one does not have to deal with singular integrals (as in the case of 
BEM) which require elaborate and time consuming algorithms for their evaluation. For this 
reason, the use of ESM in the scattering problem results in a simple and efficient code. 


FF, MHD, ATS DMN- AIAA 2009.nb | 9 


The FSC has been used to study the shielding effect of the engine noise by a wedge shaped 
airframe for which some experimental acoustic data were available [S9]. It was found that the 
agreement between the measured and predicted acoustic data was good. It is clear from this 
example that the shielding effect of fuselage and wings can be an important method of aircraft 
engine noise control. This effect can only be studied and quantified by utilizing the FSC. Note 
that the FSC provides some simple models of incident engine noise. 

The ESM requires the inversion of a large matrix whose size and memory requirement increase 
with frequency id. Typically, the memory needed for the code is proportional to to 4 and the 

computational time is proportional to co 6 . The proportionality constants are functions of the 
scattering surface area. It can be shown that on today's scientific desktop computers with one 
processor, either the size of the scattering surface must be kept small or the incident frequency 
must be assumed to be low [S8]. Fast multipole method (FMM) can be used to increase the 
upper limit of the frequency for the same scattering surface by a factor of 3 to 4 [SI 0]. However, 
even with the use of FMM, one is not able to calculate the scattering of noise from a large air- 
craft, say the size of Boeing 777, and an incident frequency of 10 khz on a scientific desktop 
computer. But the FSC has be parellelized to run on multiprocessor machines to achieve this 
goal. 

The FSC is written in FORTRAN and can be operated in batch mode. Version 2.0 is written for 
a single processor. It can also be driven by a graphical user interface (GUI) written in C++ based 
on WXWidgets. We now have two parallel versions of the FSC -v3.1 for shared memory 
platforms, and v3.2 for distributed memory machines. These versions work quite well, and will 
be delivered to NASA by the end of January 2009. They will be available to the general public 
after the invention disclosure has been filed and approved. FSC v3. 1/3.2 require very little 
computer memory. As the result of the low memory usage, practically any scattering surface 
size/frequency of interest can be handled now. But the trade-off is computational time that can 
be extensive. 

The FSC was developed for the UNIX and LINUX operating systems but it can also be adapted 
for Windows and Macintosh environments. A Users Manual is available for various versions of 
the code [S8]. 

We believe that FSC is the most advanced and efficient aeroacoustic code available today for 
computation of the scattering of incident sound by aircraft surfaces. With parallelization to run 
on multiprocessor machines, it can fully meet the needs of the aircraft and engine industry for 
product design, evaluation and noise reduction. 

3 . 3 . 1 - Input Data to FSC 

1- Excitation frequency and freestream thermodynamic variables 

2- Numerical parameters for determining the appropriate grid generation and solution strategies 

3- The scattering surface geometries (wings, fuselages, nacelles) and liner admittances 

4- Local flow variables (density, speed of sound, and Mach number vector). Code options 
include: a- no flow, b- uniform flow, c- FSC generated small perturbation compressible flow, 
and d- user supplied flow 

5- Complex values of incident acoustic pressure and acoustic velocity at FSC requested locations 
in space and time. Code options for incident sound include: a- simple point monopole(s) and 
dipole(s), b- FSC generated engine noise from nacelle alone run, and c- results from external 
noise codes capable of producing FSC type inputs. 
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6- Observer position(s): Specify arbitrary observer points or request field calculations on 
spheres, cylinders, rectangular volumes surrounding the scattering surface. An option for request- 
ing take-off or flyover footprints is provided. 

3.3.1- Output Data From FSC 

Complex values of the instantaneous acoustic pressure, velocity and intensity at the desired 
observer positions or the noise footprint regions. Graphical presentation of the output is avail- 
able in version 2.0 of the code. 

4- Plans for Further Improvements of the Codes 

It can be seen from the above presentation that Langley researchers have developed the aeroa- 
coustic technology to study the open rotor noise generation and reduction. The available codes 
are able to meet much of the need of the aircraft and engine industry for product design, evalua- 
tion and noise control. We feel that there are areas that are in need of further improvements. We 
also feel that the aircraft and engine industry have specific problems of interest that may require 
further extensions of the capabilities of our codes. We are interested to learn of these problems. 
We present below some of our research and development plans which will increase the capabili- 
ties and usefulness of our codes. 


4.1- Plans for ASSPIN Code 

ASSPIN is in an advanced stage of development and we suggest that this code be used for 
studying free open rotor noise. The reason is that to use the FSC with ASSPIN output as the 
incident sound source, we must be able to calculate the gradient of the acoustic pressure on 
all the scattering surface. The most numerically efficient method of finding the gradient is by 
an analytic method. This can be done rather easily for Formulations 1 and 1A [54] but it is 
next to impossible for Formulation 3 because of its analytic complexity and the number of 
terms. However, since the ASSPIN code was developed in nineteen eighties, there is a need 
for upgrading some of its capabilities. We recommend the following improvements to the 
ASSPIN code: 

1- Improvements of the blade geometry specification- The blade geometry of the modern 
open rotors is a complex three dimensional curved surface. We require the calculation of 
various curvatures of the surface to use in Formulation 3. This means that we need a surface 
specification with smooth second derivatives. This task was relatively easy when NACA 2D 
airfoils were used by some open rotor blade designers because analytic expressions for airfoil 
coordinates were available. However, today's blade designers use advanced airfoils optimized 
by CFD codes for their performance. We need to use sophisticated three dimensional blade 
surface specification, e.g., NURBS, which will give us surface curvature information for 
noise prediction. We have made some progress in this direction but more work needs to be 
done in this area. This feature will be added as a pre-processor to ASSPIN. 
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2- Upgrading ASSPIN to compute the noise of contra-rotating open rotors with unequal 
blade numbers and geometry in each row- Currently, ASSPIN code can handle contra-rotat- 
ing open rotors with equal number of blades on each row and identical geometry. To treat the 
general case of unequal number of blades on each row, the acoustic pressure signatures of the 
interaction tones must be computed. This involves changes to the observer time step and period 
specifications in ASSPIN code itself. Furthermore, the part of the code for the computation of 
the acoustic spectrum must be modified to account for the interaction tones. This feature will 
also give us the ability of predicting the noise of a contra-rotating open rotor with two rows of 
rotors with very different geometries (blade numbers and blade designs). 

4.2- Proposed Plans for FW - H pds code 

We have not used this code for open rotor noise prediction and what we are proposing is all new 
and requires further research. However, for two reasons we feel that work on this methodology 
must proceed with urgency because of the important capabilities it will provide us. First, we can 
include the near field nonlinear propagation effects which can be important in noise prediction. 
Second, this is the only open rotor prediction method that can be linked easily to the FSC to 
study the shielding effect and noise footprint changes by scattering from fuselage and wings (see 
Subsection 4.3, Item 2). We, therefore, propose the following research plans: 

1- Removal of spurious signals from the computed acoustic pressure signature when eddies 
and wakes cross the data surface- When turbulent eddies or the wakes from open rotors cross 
the data surface of the FW - H pds code, spurious signals appear in the acoustic pressure signa- 
ture. This is a major cause of concern for the accuracy of this method. We now outline here a 
method to alleviate the spurious signals. Our discussion refers to Figure 3 below. 
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Figure 3- Sketch of the data surface for FW - H pds code and the use of equivalent sources 
to generate the null field inside the data surface at control points. Here we show one equiva- 
lent source and one control point for simplicity. 


We first note that in the application of this code, we specifically ignore the quadrupoles outside 
the data surface because specifying these sources is not practical. We only use surface sources in 
this method. When the rotor wakes cross the data surface, we have a fluctuating quadrupole field 
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in the exterior of the data surface which generate the spurious signals. We have one important 
piece of information about the quadrupole field from the surface data. The sum of the contribu- 
tions of the quadrupoles and the surface sources inside the data surface is identically zero (null 
field) for all time. Therefore, we can use some discrete sources (generally monopoles that we 
call equivalent sources ) with unknown strengths in specified exterior or interior locations near 
the area of wake crossing to produce null field at specified control points inside the data surface. 
We must have, at a minimum, as many control points as the number of equivalent sources to be 
able to find the unknown strengths of the equivalent sources from a system of linear equations. 
Once the strengths of equivalent sources are found, the noise at any observer position in the 
exterior of the data surface is found from the surface sources as well as the equivalent sources 
and the spurious signals should disappear. 

The problem of finding the source strengths is almost certainly ill-posed of the common variety 
seen, for instance, in solving the scattering problem using equivalent sources (as in FSC). The 
handling of such ill-posed problems is understood quite well and we do not expect undue difficul- 
ties to find the strengths of the equivalent sources. 

2 - Development of an algorithm to transfer CFD data to the data surface for acoustic calcu- 
lations- This is an important recurring problem in aeroacoustics. We ask for high resolution 
CFD data for noise calculations and invariably the CFD grids do not coincide with the grids on 
the data surface for acoustic calculations. Using a surface generated by the CFD grids for acous- 
tic calculations generally does not work because we need some surface geometry information 
(outward normal, etc.) that cannot be obtained easily from such a surface. The accurate extrapola- 
tion of high resolution CFD data to the data surface presents a challenging problem that needs to 
be solved to make the FW - H pds code work satisfactorily. Such an algorithm must be imple- 
mented as a preprocessor for the code. 

4.3- Proposed Plans for the Fast Scattering Code (FSC) 

The FSC is in the process of intensive development at NASA Langley by Dunn and Tinetti. It 
has not been used for calculation of scattering of sound for an open rotor as the incident source 
of sound. However, we feel that the scattering of open rotor noise from aircraft fuselage and 
wings must be studied in the design stage for noise control. We propose the following develop- 
ments and research plans for the FSC: 

1- Search for efficient and accurate method of storage and inversion of large and dense 
matrices on multiprocessor machines- For multiprocessor machines, the selection of an effi- 
cient storage and the inversion method of large and dense matrices that ESM requires is a chal- 
lenging problem. The accuracy of the results and the execution time are both dependent on the 
method of storage and inversion of the matrix. One criterion of the search for an algorithm is 
reaching a better balance between memory and time for the very large cases. This work will be 
performed by Dunn and Tinetti for NASA Langley. 
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2- Development of a method for finding the pressure gradient of the incident sound at all 
points of scattering surface- To satisfy the boundary condition on the scattering surface 
(equality of the normal velocities of incident and scattered acoustic pressures), one needs to 
know the gradient of the incident sound. This is a very difficult problem. In the case of using 
ASSPIN for open rotor noise prediction, we have not been able to solve this problem yet. But for 
the case of FW - H pds code used to simulate the incident sound for open rotors, we have been 
able to derive two analytic results for the acoustic pressure gradient based on Formulations 1 and 
1A [54], These are known as Formulations G1 and G1A with similar characteristics as Formula- 
tions 1 and 1 A. Either result should be coded as a post- processor of the FW - H pds code and thus 
provide the input data to the FSC. 


5- Concluding Remarks 

In this paper, we started with a brief review of the contributions of Langley researchers in the 
prediction of rotating blade machinery noise. We have then presented a discussion of the open 
rotor noise prediction codes developed at NASA Langley. Three codes ASSPIN, LW - H pds and 
the Past Scattering Code (PSC) are discussed together with their capabilities. The required input 
data are pointed out and the output data are described for these codes. Pinally, some plans for the 
improvements of these codes are presented. In particular, a method based on equivalent sources 
is outlined to get rid of spurious signals in the FW - H pds code. 
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